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Abstract 

A recently formulated framework of highly-anisotropic and strongly-dissipative hy- 
drodynamics (ADHYDRO) is used to describe the evolution of matter created 
in ultra-relativistic heavy-ion collisions. New developments of the model contain: 
the inclusion of asymmetric transverse expansion (combined with the longitudinal 
boost-invariant flow) and comparisons of the model results with the RHIC data, 
which have become possible after coupling of ADHYDRO with THERMINATOR. 
Various soft-hadronic observables (the transverse-momentum spectra, the elliptic 
flow coefficient t>2, and the HBT radii) are calculated for different initial condi- 
tions characterized by the value of the initial pressure asymmetry. We find that as 
long as the initial energy density profile is unchanged the calculated observables 
remain practically the same. This result indicates the insensitivity of the analyzed 
observables to the initial anisotropy of pressure and suggests that the complete 
thermalization of the system may be delayed to easily acceptable times of about 
1 fm/c. 

PACS: 25.75.-q, 25.75.Ld, 24.10.Nz, 25.75.Nq 



1 Introduction 



The soft-hadronic observables measured in ultra-relativistic heavy-ion colli- 
sions at RHIC (transverse- momentum spectra, the elliptic flow coefficient t>2, 
and the HBT radii) are described very well by the perfect-fluid hydrodynamics 
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or by viscous hydrodynamics with a small viscosity to entropy ratio (for recent 
reviews see, for example, [Tj and [2], respectively). The implicit in this success 
is the use of a rather early thermalization time that is quite often well below 
1 fm0 The early starting time of hydrodynamics, that is identified with 
the thermalization time, is justified in a natural way if the quark-gluon plasma 
is a strongly interacting system [5|6] . 

On the other hand, there are attempts to describe the RHIC data with the 
help of models which do not assume very fast thermalization. The success 
of such models in reproducing the data may indicate that the quark-gluon 
plasma is not necessarily a strongly interacting system. An example of such 
a model has been introduced in Ref. [7], where the authors assume that the 
perfect-fluid stage is preceded by free streaming of partons, see also [8f9fTU] . 
Another example is a model discussed in Ref. [TT], where the initial stage 
consists of partons with thermalized transverse degrees of freedom only [12J. 

Recently, we have formulated a framework of highly-anisotropic and strongly- 
dissipative hydrodynamics (AD HYDRO) [T3fl4|ll5] . This model interpolates 
between a highly-anisotropic initial state (where the longitudinal and trans- 
verse pressures may be substantially different from each other) and the regime 
described by the perfect-fluid hydrodynamics. It may be used to analyze the 
effects of early anisotropic pressure on the subsequent evolution of matter. 
ADHYDRO may be treated as an important generalization of the models in- 
troduced in Refs. [7|TT] because the change from the early highly-anisotropic 
system to the subsequent locally isotropic stage is described as a continuous 
process. We note that in [TfTT] a similar change was described with the help 
of the Landau matching conditions which demand that the energy flux in the 
transition is conserved but do not demand the continuity of all components 
of the energy-momentum tensor. 

Our formulation of the model in [T5]m] has been followed by the papers by 
Martinez and Strickland [T6][TT] . where a similar approach has been derived 
from the kinetic theory with the collision term treated in the relaxation time 
approximation. In the case of purely longitudinal expansion one may show that 
the approaches introduced in [T5]rH] and [ToTfTT] lead to the same structure 
of equations that determine the evolution of pressure asymmetry and entropy 
density. The only difference may be found in the form of the entropy source. 
Interestingly, if the pressure asymmetry is small, both approaches are strictly 
equivalent and they agree with the second-order Israel-Stewart theory. In Ref. 
[T5] we have shown additionally, that the entropy source terms used in [T5]m] 
and [TfifTT] lead to practically the same results even if the initial pressure 
asymmetry is very large. 



We use the natural system of units where c = ks = h = 1. 
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In this paper, besides the boost-invariant longitudinal motion, the transverse 
expansion is also taken into account. In such a general case, the explicit con- 
nections between ADHYDRO, the underlying kinetic theory, and the Israel- 
Stewart theory have not been established yet. However, in the initial stage 
where the effects of pressure asymmetry are most significant, the longitudinal 
motion dominates the dynamics of the system, and we deal with a situation 
that is quite similar to that described in Refs. [T3fT6] . On the other hand, at 
the stage where the transverse flow becomes substantial, the pressure asym- 
metries are eliminated and the system's expansion is well described by the 
perfect-fluid hydrodynamics. Therefore, we consider ADHYDRO as a good 
candidate for the model which extrapolates well between an initial highly- 
anisotropic state and a later perfect-fluid phase o 

In this paper we use ADHYDRO and show that as long as the initial energy 
density profile is unchanged the calculated observables remain practically the 
same. This result indicates the insensitivity of the analyzed observables to the 
initial anisotropy and suggests that the complete thermalization of the system 
may be delayed to easily acceptable times of about 1 fm. In this way we find 
further support to ideas presented in [7]ITT] . 

The fact that different initial pictures lead to the same results for hadronic ob- 
servables means that the freeze-out conditions for different physical scenarios 
are similar. This suggests the formation of universal flow in heavy-ion colli- 
sions, as pointed out first in Refs. [T%][T9] . On the other hand, different initial 
conditions may affect other observables, that have not been studied here. In 
particular, the early pressure anisotropy may have an important impact on the 
directed flow, as shown recently in [20]. The study of this problem requires, 
however, a generalization of the ADHYDRO code to 3+1 dimensions, which 
is a work in progress. 



2 Definition of the model 

The ADHYDRO model is based on the following form of the energy-momentum 
tensor [T5] , 

T» v = (e + P±) UW - P ± g» u - (P ± - P^VV. (1) 



2 A natural connection between ADHYDRO and dissipative hydrodynamics may 
be achieved if the model is constructed explicitly in such a way that it interpolates 
between an initial highly-anisotropic state and a later dissipative stage described 
by the Israel-Stewart theory. This work is in progress. 
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Here s, P±, and P are the energy density, transverse pressure, and longitu- 
dinal pressure, respectively. In the case where the two pressures are equal, 
Pi — P\\ — P, we recover the energy-momentum tensor of the perfect fluid. 
The four-vector U^ 1 in flTJ is the hydrodynamic flow, 



U* 1 = >y(l,v x ,v y ,v z ), 7 = (l-w) 



2X-1/2 



(2) 



and V 1 defines the longitudinal (beam) direction, 
^ = 7,(^,0,0,1), lz = (l-v 2 z )- 1 ' 2 . 



(3) 



The four-vectors If 1 and satisfy the following normalization conditions: 



-1, u-v = o. 



(4) 



In the local- rest-frame (LRF) of the fluid element we have £/ M = (1, 0, 0, 0) 
and V* = (0, 0, 0, 1), and the energy-momentum tensor is reduced to a simple 
form, 



LRF 



( e ^ 
P ± 
Pj_ 

o p 



(5) 



Thus, the formula ([T]) allows for different pressures in the longitudinal and 
transverse directions. 



Besides the energy-momentum tensor ([T]), we introduce the entropy flux 



a' 



a 



(6) 



where a is the entropy density. We assume that e and a are functions of P± and 
P||. For massless partons the condition = gives e = 2P± + Py, however, 
in our actual calculations we use more involved expressions introduced below 
in Sect. S 

The space-time evolution of the system is governed by the equations expressing 
the energy-momentum conservation and the entropy growth, 



<Vr^ = o, 



E. 



(7) 
(8) 
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The function E represents the entropy source. The form of £ must be treated as 
an assumption that defines the dynamics of the anisotropic fluid. It is natural 
to assume that £ > 0, with £ = for P^ = P». In this way, in the case where 
the two pressures are equal, the structure of the perfect-fluid hydrodynamics 
is recovered. 

Treating E as a function of P±_ and P\\, Eqs. fl7|) and (JSJ become a closed 
system of 5 equations for 5 unknown functions: three components of the fluid 
velocity, P±, and Pu. The projections of Eq. <£7§ on U u and V u yield 



Ve = -(e + P ± )A + (P ± - P\\)U V V^V V , (9) 
Wd^ = (e + Pa_)V v VU v + (P ± - (10) 

where we have introduced the operators T> = C/^9„ and A = dJJ^. 



3 Implementation of boost-invariance 



In this paper we consider boost-invariant systems which expand arbitrarily in 
the transverse plane. Therefore, we introduce the following parameterizations 
of the four- vectors U^ 1 and V^, 



£/ M = (uo cosh?7, u x , u y , uq sinh^), (11) 



V = (sinhr/, 0,0, cosh rj), (12) 

where w , u X: and %iy circ components of the flow in the plane z = 0, satisfying 
the normalization condition 

ul-ul-u 2 y = l, (13) 



and 77 is the space-time rapidity 



, = im^. (14) 

Equations (TTTT) and ( TT2|) satisfy automatically the normalization conditions 
(j3]). In the boost-invariant case the derivative X> and the divergence A have 
the form 
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V = U^ = u d T + u ± -V ± , 

A = dJJ" = d T u + — + V ± ■ ul, 

T 



(15) 
(16) 



where uj_ = (u x ,u y ), Vj_ = (d x ,d y ), and r is the longitudinal proper time 



T 



Vt 2 - Z 2 . 



(17) 



Moreover, for the boost-invariant flow we find 



T 



Equations ( 1181) may be used to check that Eq. (TTDjl is automatically fulfilled 
for a boost-invariant system. Thus, from the energy-momentum conservation 
we obtain only three independent equations, the formula (Q which we repeat 
here for completeness 



Ve = -(e + P ± )A + (Pj_ - Py)^^^, 



(19) 



and two equations describing transverse dynamics. The latter can be chosen as 
the linear combinations: U^T^ 1 + U 2 d fl T^ 2 = and U 2 d fJ T^ 1 - Uxd^T^ 2 = 0, 
which leads us to the two expressions 



u± 



e + P ± 



V 



1 



u 2 y (e + P± 



U ± T 

(U X 8 y ~ Uyd X ) P±. 



(20) 
(21) 



Together with the entropy production equation, see Eq. (EJ), that can be put 
in the form 



Va + crA = S, 



(22) 



we have 4 equations that should determine the dynamics of 2 thermodynamics- 
like parameters (for example P± and P\\) and 2 components of the transverse 
flow (u x and u y ). 
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4 Generalized equation of state 



The equations of perfect-fluid hydrodynamics form a closed system of equa- 
tions if they are supplemented with the equation of state that specifies thermo- 
dynamic quantities as functions of one parameter, e.g., the system's tempera- 
ture (for systems with non vanishing baryon number two parameters are neces- 
sary). In the similar way, our framework requires that all thermodynamics-like 
quantities can be expressed as functions of two arbitrarily chosen parameters. 
Such relations play a role of the generalized equation of state and allow us to 
close the system of dynamic equations formulated in the previous sections. 



In our previous studies [T3~f2~l~] we have shown that it is better to switch from 
P± and Pi I to two other thermodynamics-like parameters: the non-equilibrium 
entropy density a and the anisotropy parameter x. In particular, if we deal 
with partons described by the anisotropic phase-space distribution function 
obtained by squeezing/stretching the longitudinal momentum in the Boltz- 
mann distribution by the factor y/x, we obtain 



' 2 \ 4 /3 

7rer\ 



4#o 



R(x), 



(23) 



/ttV\ 4/3 \R{x 



+ xR'(x) 



(24) 



4/3 



R{x) 



2xR'(x) 



(25) 



The function R(x) is defined by the formula [2T 



R(x) 



2tt 2 



1 + 



x 



arctan \Jx — 1 



y/x=T 



(26) 



and go is the degeneracy factor connected with internal quantum numbers of 
partons. The symbol R'(x) denotes the derivative of R(x) with respect to x, 
for x = 1 we have R'(l) = and the two pressures are equal. 

We stress that Eqs. ( !23|) - (|25|) describe the non-equilibrium state defined by the 
values of the non-equilibrium entropy density a and the anisotropy parameter 



7 



x. However, since in the equilibrium we have: E\& = 3goT 4 /TT 2 , P ic i = goT 4 /ir 2 , 
and cQj. 
form LH, 



and (Jjd = 4g Q T 3 /ir 2 [lj, Eqs. ( 123|) - (!25|) may be written equivalently in the 



£ = £ id (cr)r(x), (27) 
P± = Pid(v) [r(x) + 3xr'(x)\ , (28) 

P y = P id ( a ) [ r (x) - 6zr'(x)] , (29) 

where £id(c) and Pid(o") are equilibrium expressions for the energy density 
and pressure, and r(x) = 7r 2 R(x)/(3go). One can easily notice that for x = 1 
we have r(l) = 1 and r'(l) = 0, hence Eqs. (l27]) - (l29l are reduced to the 
equilibrium expressions. 

As stated above, the structure of Eqs. (j23])-([25]) and Eqs. (E7J) ([29D follows 
from the special form of the momentum distribution. This form is very much 
likely to be realized at the very early stages of the collisions. On the other 
hand, at later stages of the collisions, the produced matter approaches local 
equilibrium described by the appropriate equation of state. In other words, as 
the system becomes more isotropic, i.e., as the parameter x tends to unity, 
the local properties of the fluid become closer to those characterized by the 
realistic equation of state. This argument convinced us to use, in the numerical 
calculations presented below, the following expressions 

e = e qgp (a)r(x), (30) 
P± = P qgP (cr) [r(x) + 3xr'{x)} , (31) 

P\\=P^)Hx)-Qxr'(x)], (32) 

where e qgp (a) and P qgp (a) characterize the realistic equation of state for van- 
ishing baryon chemical potential, as it has been constructed in Ref. [22J, and 
the function r(x) is the same as that appearing in Eqs. (!27j) - (l29l) . 

It is necessary to admit that Eqs. (13"U|) - (!3"2"|) have no direct microscopic expla- 
nation. The main motivation for this form is that it has two attractive limits. 
At the very early stage the system consists most likely of massless partons 
and Eqs. fl30l) - fl32l) approach Eqs. ff2~T|) - fl2"9"]) . since the realistic equation of 



3 For simplicity, Eqs. (I23p - (l25p have been introduced for the squeezed/stretched 
classical (Boltzmann) distribution. A generalization to the squeezed/stretched quan- 
tum (Bose-Einstein or Fermi-Dirac) distributions has been given recently in [15j . 
This generalization introduces an irrelevant factor. 
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state approaches the Stefan-Boltzmann limit (although deviations from this 
limit are noticeable). On the other hand, when the system becomes isotropic, 
x = 1 and Eqs. (l3"0"]) - ([3"2"]) are reduced to the equation of state used in standard 
hydrodynamics. 



5 Entropy source 



In this paper, following our original formulation of the model in Ref. 
use the following form of the entropy source 

. o-v^ (33) 



X T e q 



The quantity r eq is a timescale parameter. The form ( |33|) guarantees that 
E > and E(er, x = 1) = 0. More arguments for such a particular form 
of E are given in [13]. They connect x with the ratio of the transverse and 
longitudinal temperatures. 

In the case of purely longitudinal expansion and for small deviations from 
equilibrium one can show that Eq. ( 133]) leads to quadratic dependence of the 
entropy source on the variable £ = 1 — x and this dependence is compati- 
ble with the Israel-Stewart theory (where the entropy production depends on 
the viscous stress squared) and with the Martinez-Strickland model [16J. For 
more general situations, the equation ( |33l) should be treated as one of the 
assumptions defining our model. Other forms of the entropy source have been 
analyzed in [15]. All those forms lead to similar numerical results. 

We emphasize that the structure of the entropy source is an external input 
for the anisotropic hydrodynamics. Especially, in the region where the asym- 
metries are large and no correspondence to dissipative hydrodynamics can be 
found. In this context, it is interesting to obtain any hints about E for large 
x from the microscopic models of particle production or from the AdS/CFT 
correspondence. 

In our numerical calculations we adopt the value r cq = 0.25 fm. This choice is 
suggested by our previous calculations [T3~|14|I15] . where we have shown that 
it leads to almost complete thermalization of matter at the proper time of 
about 1 fm (for highly anisotropic initial conditions studied here). Our earlier 
calculations have shown that reaching the perfect-fluid regime at r ~ 1 fm is 
necessary to obtain a good agreement with the data 0- 

4 If the free-streaming or transverse-hydrodynamics stage is longer than 1 fm, it is 
difficult to achieve a good description of data, as discussed in Refs. |7|llj . 
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The value r eq = 0.25 fm corresponds to r™ 11 = (2/15)r eq ~ 0.03 fm used in the 
collision term by Martinez and Strickland [T5]fT6] . Such a small value of Tg° u 
indicates that the plasma is indeed a strongly interacting system. However, our 
example shows that the average collision time in a highly-anisotropic plasma 
and the plasma isotropization time may be quite different. 



6 Anisotropy evolution 



If both the generalized equation of state and the entropy production term 
are defined, one can derive a compact expression for the time evolution of 
the anisotropy parameter x. Of course, this expression should be considered 
together with other dynamic equations. Nevertheless, its form turns out to be 
useful in the analysis of the time evolution of the system. 



In the case described by Eqs. fl30|) -f l32|) . we find 



Vx = ^ - A) - ( 1 + (34) 

£qgp V T / \ £qgp / r eq 



where 



= M<W£)!. (35 ) 
r[x) 



If 3P qgp = £ qgp , which is the case realized by Eqs. (|27I) - ([25I) . Eq. (134"|) is reduced 
to the form 



Vx = x (— -A]- —H(x) . (36) 



\ r ) 3r, 



ex. 



For the purely longitudinal boost-invariant motion, uq = 1 and A = 1/r, and 
Eq. ( 136]) simplifies to an ordinary differential equation for x, 

dx = 2x_ AH(x) 
dr r 3r eq 



We shall come back to the discussion of Eq. (1371) in Sect. 18.11 
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7 Initial conditions 



The initial conditions for the evolution are defined by four functions: cr(r , xj_), 
x ( T o, x _l)> u x (r , xj_), and u y (T , xj_), where r is the initial proper time. In the 
numerical calculations we assume r = 0.25 fm. We also assume that there 
is no transverse flow present initially, therefore we set u x (t , xj_) = and 
%(to,x_l) = 0. 

For the initial anisotropy x(ro,xj_), which we take as independent of xj_ and 
denote simply as Xq, we consider three different options: xq = 100, xq = 1, 
and Xq = 0.032. The case x — 1 is, of course, the closest to that described by 
standard perfect-fluid hydrodynamics. 

The case Xq = 100 corresponds to the initial situation where the transverse 
pressure is much larger than the longitudinal pressure. In this case the momen- 
tum shape is oblate; the momentum distribution is stretched in the transverse 
direction and squeezed in the longitudinal direction to the beam. This type 
of the initial conditions is considered in the Color Glass Condensate (CGC) 
approach where the distribution functions in the longitudinal direction are 
described by the Dirac delta function, 8(p\\) at z = [2"3"f2"l] . see also Ref. [2"5] . 
The values xq < 1 correspond to the prolate momentum shape. This type of 
the initial conditions have been analyzed, for example, in Refs. [26|27] . 

The value xq = 0.032 has been chosen, since r(100) ~ r(0.032), and for fixed 
energy density the entropy densities in the cases x = 100 and x = 0.032 are 
the same. Therefore, the case x = 0.032 may be considered as a counterpart 
of the case x$ = 100 with inversed role of pressures. 

In all considered cases we assume that the initial energy density in the trans- 
verse plane, Eo(to,x±), is proportional to the normalized density of sources 
p(b, Xj_), where b is the impact parameter corresponding to a given centrality 
class, 

£o(t ,x_l) = £ip(6,x_L). (38) 



The parameter e- x is the initial energy density at the center of matter in most 
central collisions. Its value may be estimated from the standard hydrody- 
namic calculations [1J. We use the value e\ = 86.76 GeV/fm 3 . If the consid- 
ered matter was in equilibrium, its temperature would be equal to 485 MeV. 
The normalized density of sources p(b, xj_), is constructed as a combination 
of the wounded- nucleon density pwip, xj_) and the density of binary collisions 
PbQ>,x±) [28], 

1 - a 

p(b, x ± ) = — — p w (b, x ± ) + ap B (b, x ± ) , (39) 



11 



(40) 



The distributions pw (b, xj_) and ps (b, xj_) are calculated for a given centrality 
class from the Glauber model in the optical approximation. Following the 
PHOBOS studies of the centrality dependence of the hadron production [2S] 
we take a = 0.14. In the limit a — > our assumption (|39|) is reduced to the 
wounded nucleon model [30]. 



By fixing both the initial energy density and the initial anisotropy parameter, 
we determine the initial entropy density profile from Eq. ( 130]) . namely 



0-(t o , x ±) = £ 



-l 

;qp 



£j/j(6,X ± ) 

r(x ) 



(41) 



Here e~^(e) is the inverse function to e gqp (a). The 6-dependence displayed 
on the right- hand- side of Eq. (I41jl induces centrality dependence of the initial 
entropy density profiles. 



8 Results 



8.1 Central collisions 



Let us first discuss our results obtained for the centrality class c = — 5%, 
which corresponds to the impact parameter b = 2.26 fm. In our all calculations 
we use the value e\ = 86.76 GeV/fm 3 . This implies that the initial entropy den- 
sity at the center equals <t ; = cr(r ,x^ = 0) = 227.6 fm -3 for the case x = 1, 
and o\ = 151.7 fm -3 for the cases Xq = 100 and Xq = 0.032. The same initial 
energy density, e\ = 86.76 GeV/fm 3 , is assumed in the perfect-fluid calculation 
based on the LHYQUID code [3~Tf3"2] . which serves as a reference point. Since 
the perfect-fluid case is obtained formally by setting x = in ADHYDRO, the 
initial entropy density in the perfect-fluid calculation equals a- x = 151.7 fm -3 
and the initial temperature is 7] = 485 MeV. 

Figure [1] shows the time dependence of the entropy density at the center of the 
system, i.e., the function o~(t, xj_ = 0), in the time interval 0.25 fm < r < 2 fm, 
for three different initial values of the anisotropy parameter: xq = 100 (dashed- 
dotted line), xq = 1 (dotted line), and xq = 0.032 (dashed line). The entropy 
density increases significantly at the initial stage of collisions in the two cases 
exhibiting strong initial asymmetry, i.e., for xq = 100 and xq = 0.032. After 
reaching the maximum at r = 0.3 fm, <r(r, xj_ = 0) starts to decrease. For 
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300 




0.30 0.50 0.70 1.00 1.50 2.00 

r [fm] 

Fig. 1. Time dependence of the entropy density at the center of the system, 
<t(t, = 0), for three different initial values of the anisotropy parameter: xq = 100 
(dashed-dotted line), xq = 1 (dotted line), and xq = 0.032 (dashed line). The results 
have been obtained for the centrality class c = — 5%. The corresponding result of 
the perfect-fluid hydrodynamics is denoted by the solid line. 

r > 0.5 fm, the entropy density scales approximately as 1/r. This behavior 
reflects the approximate perfect-fluid behavior reached for significantly large 
times GO- In the case Xq = 1, the entropy density decreases in the very much 
similar way as in the perfect-fluid case denoted by the solid line. Interest- 
ingly, in all considered cases the final entropy densities are very much similar. 
This result suggests that the final particle multiplicities may be practically 
the same. 

The behavior of the entropy density may be discussed from a complementary 
point of view if we analyze the function S(r, xj_ = 0). The time dependence of 
the entropy source is shown in Fig. [2j The notation is the same as in Fig. [TJ In 
the case x = 100 (dashed-dotted line) the initial entropy production is large. 
At later times E(r, xj_ = 0) becomes smaller, since the pressure asymmetry 
decreases; see Fig. [3] where the corresponding time dependence of x(r, x^ = 0) 
has been shown. In the case xq = 0.032 (dashed line) the initial entropy 
production is also large but it drops to zero at the time when x passes unity. 
Later it increases and becomes similar to that found for the case Xq = 100. If 
the initial asymmetry is 1 (dotted line), the entropy production is generally 
small because x stays close to unity during the whole evolution. 

We stress that the results obtained with ADHYDRO where x has been initially 
set equal to zero are slightly different from the results obtained within perfect- 
fluid hydrodynamics. The origin of this difference may be traced back to the 

5 The scaling a oc 1/r is strictly valid for the longitudinal, boost-invariant expan- 
sion. The transverse expansion induces corrections to this behavior of about 10%. 
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Fig. 2. Time dependence of the entropy production at the center of the system, 
^( T ) x _i_ = 0))) f° r three different initial values of the anisotropy parameter, notation 
the same as in Fig. [T] (c = — 5%). 
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Fig. 3. Time dependence of the anisotropy parameter at the center of the sys- 
tem, x(t, xj_ = 0), for three different values of xq, notation the same as in Fig. [1] 
( c = 0-5%). 

terms in the evolution equation for x, which describe free streaming. In the 1+1 
boost-invariant case with Eqs. (!27|) - (l29|) . which is easy to analyze explicitly 
[13], the time dependence of x is determined by the equation (157|) . 

The first term on the right-hand-side of (|3ZD describes the effects of free stream- 
ing, while the second term is responsible for thermalization/isotropization. For 
small initial values of x, the first term dominates the very early dynamics and 
x starts to grow with time. However, as the evolution time increases, the first 
term becomes negligible, the second term takes over and determines an ap- 
proach to equilibrium. For large initial values of x the second term dominates 
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the dynamics of matter from the very beginning. If we assume x = 1 at the 
beginning of the evolution, the second term vanishes, since H(l) = 0, and x 
increases. This effect is important only if r is small, which agrees with the idea 
that the free streaming of partons may be relevant only at the beginning of 
the evolution of matter, when the system is very much dense and the effects 
connected with asymptotic freedom may play a role. 

Having solved numerically the ADHYDRO equations for given initial con- 
ditions we determine the freeze-out hyper surf aces. We follow the same pro- 
cedure as in standard 2+1 hydrodynamics. The only difference is that we 
exclude times smaller than 1 fm, where the system is substantially anisotropic 
and consists of gluons rather than hadrons that can be freely emitted. For 
t > 1 fm, the matter is practically in local equilibrium and we may specify 
the freeze-out condition by fixing the freeze-out temperature Tf. In our calcu- 
lations we assume Tf = 150 MeV, the value which turned out be successful in 
the description of the RHIC data in our earlier approaches. If the freeze-out 
hypersurface is determined, we use it as an input in THERMINATOR 2 [34, 35] to 
generate physical events. 

Since the analysis of hadronic abundances indicates that the baryon number 
density is not exactly zero at midrapidity in Au+Au collisions at ^/s^ = 200 
Ge V, we include small non-zero values of the chemical potentials [36] . 

In Fig. H we show the transverse- momentum spectra of pions (upper part), 
kaons (middle part), and protons (lower part). We use the same notation as 
in the previous Figures: x = 100 (dashed-dotted line), x — 1 (dotted line), 
Xq = 0.032 (dashed line), LHYQUID (solid line). One can see that all versions 
of the calculations yield very similar results for the spectra. A small persistent 
difference can be noticed — the results with xq = 100 are the highest, while 
those with x = 0.032 are the lowest. In particular, the spectra obtained with 
xq = 100 agree better with the data in the high p± region. Nevertheless, 
the results obtained with different values of x are close to each other and 
to the standard 2+1 hydrodynamic result obtained with LHYQUID. This 
means that the strong early anisotropic behavior have small impact on the 
observables. This, in turn, leads us to the conclusion that the data does not 
exclude a possibility that such an anisotropic stage can exist for about 1 fm at 
the beginning of the collision. It is also preferable, that the initial transverse 
pressure is larger than the longitudinal pressure. 

Figure E shows the HBT radii calculated with FEMTO THERMINATOR, 
a part of the THERMINATOR 2 package [33]. FEMTO THERMINATOR 
implements the methods described in more detail in Ref. [38]. We observe 
again a small dependence of the model results on the initial value of x. The 
agreement between the model results and the data is very good for R out and 
-Riong- On the other hand, the model results for -R S id e are about 10% below the 
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Fig. 4. Transverse- momentum spectra of pions (upper part), kaons (middle part), 
and protons (lower part) for Au+Au collisions at - v /snn = 200 GeV and the cen- 
trality class c = — 5%. Notation the same as in the previous Figures. The data 
are taken from |33j . 

experimental values. Clearly, the HBT puzzle is not completely eliminated in 
our approach. The presence of the initial, strongly anisotropic stage does not 
help to improve the agreement with the data that has been already obtained 
with the LHYQUID calculation. Nevertheless, the overall agreement with the 
model results and the data is quite satisfactory. 
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Fig. 5. HBT radii calculated with ADHYDRO and LHYQUID. The model results 
for Au+Au collisions at ^/snn = 200 GeV and the centrality class c = — 5% are 
compared to the RHIC experimental data |37j . 

8.2 Non-central collisions 

In this Section we analyze non-central collisions. Two centrality classes are 
considered: c = 20 — 30% and c = 20 — 40%, which correspond to the 
impact parameter b = 7.16 fm and b = 7.84 fm, respectively. With the 
value £i = 86.76 GeV/fm 3 used in Eq. (14T]) . the 6-dependence of the source 
density implies the following values of the central initial entropy density: 
G\ = 183.1 fm -3 for the case Xq = 1, and o x = 120.3 fm -3 for the cases Xq = 100 
and Xq = 0.032 if c = 20 — 30%; = 171.3 fm -3 for the case x = 1, and 
a; = 113.8 fm" 3 for the cases x = 100 and x = 0.032 if c = 20 - 40%. 

In Fig. |6] we show our model results for the transverse-momentum spectra of 
pions (upper part), kaons (middle part), and protons (lower part), obtained 
for the centrality class c = 20 — 30% and compared with the RHIC data. 
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Fig. 6. Transverse- momentum spectra of pions (upper part), kaons (middle part), 
and protons (lower part) for Au+Au collisions at - v /snn = 200 GeV and the cen- 
trality class c = 20 — 30%. Notation the same as in previous Figures. The data are 
taken from |33j . 

One can notice the overall good agreement between the model results and the 
data. Interestingly, the scaling of the initial entropy based on Eq. ( HTI) ensures 
the good normalization of the spectra. Similarly to the central collisions, the 
case with xq = 100 describes slightly better the high p± tails, although the 
differences between the model and the data in this region are larger than 
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Fig. 7. HBT radii calculated with ADHYDRO and LHYQUID. The model results 
for Au+Au collisions at y / «NN = 200 GeV and the centrality class c = 20 — 30% are 
compared to the RHIC experimental data |37j . 



in central collisions. In any case, for p± < 2 GeV, where the hydrodynamics 
model is expected to work, the agreement with the data is very good. 

Figure [7] shows the HBT radii for the centrality class c = 20 — 30%. The 
agreement with the data is similar as in the case of the central collisions. 
Again, no significant effects of the early anisotropic phase are seen. 

Finally, in Fig. [H] we show the differential elliptic flow i>2(p_l) of pions and kaons 
for the centrality class c = 20 — 40%. The results obtained with different initial 
values of the anisotropy parameter coincide with each other. They are also 
consistent with the LHYQUID result. For p± < 2 GeV the model results agree 
with the data. For larger values of transverse-momentum, the perfect-fluid 
hydrodynamics overshoots the experimental data, which is a known effect. 
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Fig. 8. The elliptic flow of pions and kaons in Au+Au collisions at ^/snn = 200 
GeV and the centrality class c = 20 — 40%. The model results are compared to the 
RHIC experimental data [39|. 

9 Conclusions 



In this paper, a recently formulated model of highiy-anisotropic and strongly- 
dissipative hydrodynamics (ADHYDRO) has been used to study the effects 
of strong initial anisotropy of pressure on soft hadronic observables stud- 
ied at RHIC (transverse- momentum spectra, the elliptic flow coefficient t>2, 
and the HBT radii). We have found that the initial conditions with different 
anisotropies lead to similar results, provided the initial energy density profile 
in the transverse plane is the same. This result agrees well with earlier find- 
ings where the perfect-fluid stage starting at about 1 fm was preceded by free 
streaming of partons or by expansion of matter thermalized only in the trans- 
verse direction. The main conclusion from those studies is that the complete 
thermalization of matter may take part at the times of about 1 fm. In order to 
reproduce the elliptic flow of pions, it is not necessary to assume that matter 
thermalizes within a fraction of one fermi, as the decrease of spatial eccen- 
tricity of the reaction zone is compensated by the building of the transverse 
flow in an early non-equilibrium stage (described as free streaming, transverse 
hydrodynamics, or highiy-anisotropic and strongly-dissipative hydrodynamics 
used in this paper). 



10 Appendix: Integral form of the conservation laws 



The conservation laws (J7J) and dHJ can be integrated over the spatial transverse 
coordinates. This leads to the global conservation laws. In the case of Eq. (|S]) 
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we obtain 

J d 2 x± cr(r, x^)-u°(r, xj_)r =T J d 2 XjJl(T, xj_). (42) 

If the entropy source term E(r, xj_) vanishes, the left-hand-side of Eq. ( 14"2"|) is 
zero. This implies that the entropy per unit rapidity is conserved. In the case 
of Eq. ([8]) we define the quantity e±, 

e±(T,xj.) = e(r,xj.)^(r,xj.) + P±(r, x ± )u 2 L (r, xj. (43) 
The integration of (JTj) yields 

| d 2 x ± [e ± (r,x ± )r] = -| rf 2 x ± P||(r, x ± ). (44) 



If the longitudinal pressure P\\(t, xj_) vanishes, no work is done in the longitu- 
dinal direction, and Eq. (jUJ) leads to the conservation of the energy per unit 
rapidity. Equations (j42p and fl44|) are useful for checking the numerical code. 
We have found that in our calculations they are satisfied with the accuracy of 
about 0.001%. 
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